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Abstract 

We study the time evolution of the ASEP on a one-dimensional torus with 
L sites, conditioned on an atypically low current up to a finite time t. For a 
certain one-parameter family of initial measures with a shock we prove that 
the shock position performs a biased random walk on the torus and that the 
measure seen from the shock position remains invariant. We compute explic- 
itly the transition rates of the random walk. For the large scale behaviour 
this result suggests that there is an atypically low current such that the op- 
timal density profile that realizes this current is a hyperbolic tangent with 
a travelling shock discontinuity. For an atypically low local current across a 
single bond of the torus we prove that a product measure with a shock at an 
arbitrary position and an antishock at the conditioned bond remains a convex 
combination of such measures at all times which implies that the antishock 
remains microscopically stable under the locally conditioned dynamics. We 
compute the coefficients of the convex combinations. 

Keywords: Asymmetric symmetric exclusion process; large deviations; shocks; 
antishocks. 
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1 Introduction 



The Asymmetric Simple Exclusion Process (ASEP) on Z [HEj describes the Marko- 
vian time evolution of identical particles on Z according to the following rules: (1) 
Each particle performs a continuous time simple random walk on with hopping rate 
c r to the right and q to the left; (2) The hopping attempt is rejected when the site 
where the particles tries to move is occupied. Without loss of generality c r > q 
is assumed. An important family of invariant measures for this process are the 
Bernoulli product measures parametrized by the particle density p. The stationary 
current j* is then j* = (c r — q)p(1 — p). For c r = q (symmetric simple exclusion 
process) the dynamics is reversible and one has j* = 0. The process may also be 
defined on a finite integer lattice with various types of boundary conditions which 
may lead to non-vanishing stationary currents even in the symmetric case. The 
presence of a stationary current indicates the lack of reversibility of the dynamics 
and therefore one is interested not only with its stationary expectation j*, but also 
its fluctuations. 

Starting with the seminal papers [31 H] the large deviation theory for the ASEP 
has been developed in considerable detail. Of particular interest in this context is the 
time-integrated current Jkif) = J£(t) — J^(t) across a bond (k, k + 1), where J£(t) 
is the number of jumps of particles from k to k + 1 up to time t and analogously 
J^(t) is the number of jumps from k + 1 to k up to time t, starting from some 
initial distribution of the particles. In a periodic system with L sites, i.e., on the 
torus Tl '■= Z/LZ, a related quantity of interest is the time-integrated total current 
J(t) = J2kei L Jk(t) which is intimately related to the rate of entropy production 
[5], [6] . We denote its distribution by Zj(t) := Prob [ J(t) = J] with J G Z. One also 
consideres the (local) mean current jkif) = Jk(t)/t and the (global) mean current 
density j(t) = J(t)/(Lt). For the Bernoulli product measure one has ]im t -H X1 jk(t) = 
lim^oo j (t) = j*. The probability to observe for a long time interval t an untypical 
mean j ^ j* is exponentially small in t. This is expressed in the large deviation 
property [7] Zj(t) oc exp (-f(j)Lt) where f(j) is the rate function which plays a 
role analogous to the free energy in equilibrium statistical mechanics. Indeed, in 
complete analogy to equilibrium one introduces a generalized fugacity y = e s with 
generalized chemical potential s and also studies the generating function Y s (t) = 
(y J ^) = J2jU" J Zj(t). The cumulant function g{s) = lim 4 _ 5 . 00 lny s (t)/(Lt) is the 
Legendre transform of the rate function for the mean current, g(s) = maxj\js—f(j)]. 
The intensive variable s is thus conjugate to the mean current density j. 

Recently a focus of attention has been on the spatio-temporal structure of the 
process conditioned on realizing a prolonged untypical behaviour of the current. A 
convenient approach to study this rare and extreme behaviour is to consider the 
process in terms of the conjugate variable s. Fixing some s ^ corresponds to 
studying realizations of the process where the current fluctuates around some non- 
typical mean. We shall refer to this approach as grandcanonical conditioning, as 
opposed to a canonical condition where the current J(t) would be conditioned to 
have some fixed value J. For the weakly asymmetric exclusion process (WASEP) 
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on Tl where c r — q = v/L is small, it turns out in the hydro dynamic limit that 
for any strictly positive s (i.e., for any current j > j*) the optimal macroscopic 
density profile p(x, t) that realizes such a deviation is time-independent and flat [5]. 
Thus on macroscopic scale the conditioned density profile is equal to the typical 
(unconditioned) density profile. However, a recent microscopic approach using a 
form of Doob's /i-transform has revealed O [10] that for large s the ASEP with arbi- 
trary strength of the asymmetry c r — q exhibits interesting stationary correlations. 
Even more remarkably, the process undergoes a phase transition from the typical 
dynamics in the KPZ universality class with dynamical exponent z = 3/2 to ballistic 
dynamics with z = 1. 

For currents below the typical value (s < 0) the large deviation approach for the 
WASEP yields a different phenomenon [8]: For currents close to the typical current 
(but j < j*) the optimal macroscopic density profile p(x,t) is time-independent and 
flat as is the case for any s > 0. However, below some critical value s c (i.e. below 
some critical j c < j*) the flat profile becomes unstable and a travelling wave of 
the form p(x — vt) develops. In the limit v — > oo (expected to correspond to finite 
asymmetry in the ASEP) the optimal profile in this regime is predicted to have 
the form of a step function with two constant values p\ and pi- This is a profile 
consisting of a shock discontinuity where the density jumps from pi to pi > p\ 
at a position x\(t) and an antishock where the density jumps from p2 to p\ at 
position X2{t). Hence, conditioning on a sufficiently large negative deviation of 
the current j < j c induces a phase separation in the WASEP into a low density 
segment of length r = X\ — x<i with density p\ and a high density segment of length 
L — r = x 2 ~ X\ with density p 2 (positions and distances taken modulo L). A similar 
phase separated optimal profile, but without drift, has recently been obtained for 
atypically low activities K{t) = J + (t) + J~{t) in the SSEP [TL]. 

In order to understand this phenomenon better we develop in this paper a mi- 
croscopic approach for the ASEP for atypically low currents, but without restriction 
to the case of weak asymmetry. For a periodic system, as discussed above, one may 
expect in the long-time regime a travelling wave. In order to get insight into the 
microscopic structure of the travelling wave we focus on short times and address the 
question how for global conditioning a certain initial distribution which already has 
a shock evolves at some later finite time t. We also consider local conditioning of 
the current on a single bond. For this case we study the evolution for finite times t 
of a shock/ antishock pair with densities p\^. Our rigorous analysis is an adaptation 
of the algebraic techniques developed in [12] and extended recently by Imamura and 
Sasamoto to study current moments of the ASEP [13]. This approach - which has 
its probabilistic origin in the self-duality of the ASEP [HJ [T5j [16] - constrains our 
approach to a specific family of initial shock measures which are certain functions 
of the hopping asymmetry. However, going from the weakly asymmetric case all the 
way to the totally asymmetric limit allows us to cover a wide range of scenarios. 

The paper is organized as follows. In Sec. 2 we describe informally the well- 
established, but not so widely known tools required for studying grandcanonically 
conditioned dynamics. In Sec. 3 we introduce a new family of shock measures which 



3 



are relevant for the grandcanonically conditioned time evolution of the ASEP. In 
Sec. 4 we state and prove the main result and in Sec. 5 we conclude with some 
further comments on a generalization of the present approach to more general shock 
measures and on the macroscopic large deviation results of [8]. 

2 ASEP conditioned on a low current 

We consider the ASEP on the one-dimensional torus with L sites. For conve- 
nience we take L even and denote sites (defined modulo L) by an integer in the 
interval — M + 1, M where M := L/2. We shall identify jumps to the "right" with 
clockwise jumps. We introduce the hopping asymmetry 



and parametrize the densities p by the fugacity 



The shock/ antishock initial distributions considered below are characterized by two 
fugacities Z\p satisfying 

Condition S: 

7 = " 2 < 3 » 

By convention we shall assume 1 < q < oo, i.e. we consider preferred hopping to 
the "right" (clockwise) and we exclude the symmetric case q — 1. 

Grandcanonically conditioned Markov processes may be studied in the spirit 
of Doob's /i-transform [IT], as was done for the ASEP conditioned on very large 
currents in [13 [9], [10] , see below for the precise definition. To this end, we follow our 
earlier work [12] and we employ a matrix formulation to define the process. This 
formulation of Markovian dynamics, formalized in probabilistic terms for exclusion 
processes in [TJ5] and described informally in detail in [2], allows for an alternative 
convenient formulation particularly of conditioned interacting particle systems, as 
discussed in [6]. For self-containedness we briefly explain the approach and introduce 
the necessary notation. 



In the context of interacting particle systems with state space V the Markov genera- 
tor L acting on cylinder functions f{rj) of the configuration r\ G V is usually defined 
through the relation 




(1) 



2.1 Matrix formulation of the ASEP 




(4) 
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where uy^ is the transition rate from a configuration 77 to 7/. In particular for a 
probability measure fj,(t) we have 



f t (f),= (Lf), (5) 

where ( • ) denotes expectation w.r.t. /x In particular, taking / to be the indicator 
function l v on a fixed configuration 77 yields the master equation 

^fa; t ) = Yl K^'MV; - 1)\ (6) 

77' GV 

for the time evolution of the probability t) of finding the configuration 77 at 
time t. A stationary distribution, i.e., an invariant measure satisfying dfi/dt = is 
denoted by //* and a stationary probability of a configuration 77 is denoted fJ,*(r]). 

A convenient way to write the master equation ([6]) in a matrix form is provided 
by the so-called quantum Hamiltonian formalism [191 E]- The idea is to assign to 
each of the possible configurations 77 a column vector 1 77 ) which together with the 
transposed vectors ( 77 1 form an orthogonal basis of a complex vector space with 
inner product ( rj \ ) = 5 ViV i . Here S VtV i is the Kronecker symbol which is equal to 
1 if the two arguments are equal and zero otherwise. Therefore a measure can be 
written probability vector 

\Kt)) = ^Kv,t)\v)- (7) 

whose components are the probabilities n(rj;t) = (r)\fi(t)). The bra-ket notation 
for vectors is an elegant tool borrowed from quantum mechanics. 

Using a standard argument (see for example Chapt. XVII of Feller [20J), for all 
t > 0, the master equation ([6]) then takes the form of a Schrodinger equation in 
imaginary time, 

j t \»(t)) = -H\Kt)) (8) 

with the formal solution 

|M*)>=e- m |/i(0)> (9) 

reflecting the semi-group property. The off-diagonal matrix elements H ViV t of the 
matrix H are the (negative) transition rates w V:V ' and the diagonal entries H v>v are 
the sum of all outgoing transition rates w v ' jV from configuration 77. 
For the study of expectations we also define the row vector ( s | 

( S |:=(l,l,...,l)= Yl (l\ ( 10 ) 

r?£{0,l} L 

which we call the summation vector. By conservation of probability we have for any 
probability vector (s\/i(t)) = J2 v ev M 7 /! = 1- This implies (s\H — 0, i.e., the 
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row vector ( s | is a left eigenvector of H with eigenvalue 0. A right eigenvector with 
eigenvalue is an invariant measure | fi* ) . 

Expectation values (f) = ^ f{jf)^{ji) of a function f{rf) are obtained by taking 
the scalar product (s\f\ fi) of the diagonal matrix 

r?GV 

where | rj ) ( rj | is used as a shorthand for | rj ) ® ( rj | as is standard in the quantum 
mechanics literature. Moreover, we introduce for t > the non-diagonal matrices 

fit) = e m fe~ m (12) 

and find for any initial measure /i(0) the useful identity 

(/(*)> = (8\f\rtt)) = (8\f{t)\p(p). (13) 

Correspondingly for joint expectations at different times > t,i > one has 

( fn(tn) ■ ■ ■ h(t 2 )fl(tl) ) = (S \f n (t n ) . . . / 2 (t 2 )A(ti) | MO)- (14) 

Notice that (EJ) translates into d/(dt)(/(t) } = —(fH). As there will be no danger 
of confusing L and H we shall somewhat loosely refer also to H as generator of the 
process, since we can identify H with usual generator acting on cylinder functions 
/ by its action to the left on ( s \ fH. 

For the ASEP the vector representation of the state space V = {0, 1} L is con- 
veniently done using a tensor basis. The tensor product is denoted by ®, and A® k 
denotes the fc-fold tensor product of A. The superscript T indicates transposition. 
A configuration rj G {0, 1} L will be represented by the vector |7/) 6 C 2 which is 
defined in the following manner 

1 77 ) = I r]{-M + 1))®| rj(-M + 2) ) ® ■ ■ • <g> | rj(M) ) (15) 

where for each i = —M + 1, . . . , M, | r](i) ) = (0, 1) T , if the z-th site in the configu- 
ration 7] contains a particle, and | r)(i) ) = (1,0) T , otherwise. Observe that for any 
configuration 77, the corresponding vector | rj ) has 1 at one of its components and 
at all others. 

In order to construct the matrix representation of the generator H for the ASEP 
we introduce the three 2-by-2 matrices 

a+ := ( °o I ) ' ^ := ( 1 ) ' " := ( 1 ) ' (16) 

Then, denoting by 1 the two-dimensional unit matrix, we introduce for each fcGT 
a± ^l®...®!®^®!®---®!, n fe :=l®---®l®n®l®---®l (17) 
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where a non-identity matrix appears at the (M + k)-th position, counting from the 
left to the right. We also introduce v := 1 — h and Vj~ = t® L — h\.. 

For readers unfamiliar with this construction we point out some simple proper- 
ties of these matrices and their action on vectors. Notice that (1, 0)cr + = (0, l)h = 
(0,1), (0,l)o- = (1,0), (l,0)n = (l,0)o- = (0,l)a+ = (0,0); and analogously, 
a + (0, if = (1, Of, 0-(l, Of = n(0, if = (0, if, n(l, Of = a+(l, Of = a"(0, if - 
(0,0). These identities and the representation ()15p yield the following properties 
of and n^. If a configuration 77 does not have a particle at the site k then 
I r]') : = cr^T 1 77 ) corresponds to the configuration that coincides with 77 on \ {k} 
and has a particle at k; if to the contrary 77 has a particle at k then 0^)77) = 0, 
i.e., the vector with all components equal to 0. If a configuration 77 has a particle 
at the site k then | 77* ) := cr^~ j 77 ) corresponds to the configuration that coincides 
with 77 on \ {k} and does not have a particle at fc; if to the contrary 77 does 
not have a particle at k then 0^1 77) = 0. Using the equality (crj^f = and the 
above notations, we have that if 77 does not have a particle at k then ( 77 \ = ( 77' | 
and (77 jo"^ = 0, while if 77 has a particle at k then (77 = and (77 = (77* |. 
Accordingly, a^T and a\ are called 7j/ie 'particle creation/ annihilation operators. The 
operator n^. is called the number operator, when applied to |?7), it returns I77), if 
there was a particle at the site k, and results in otherwise, i.e. we have 

am = ,mi,) = { !"> , us) 

%l>(> = (i->fW)h} = { J,"' H^jli ■ (W) 

With these definitions the generator of the ASEP becomes 

M 

H := - ~ Wi+i) + ^0* r °i+i - «<Wi+i)] • (2°) 

j=-M+l 

For an arbitrary distribution /i(0) on {0, 1} L and its vector representation | /i(0) ) 
as defined above we denote by | /i(/j) ) the vector representation of the distribution 
of ASEP at time t, starting from /i(0). We point out that in the tensor basis a 
product measure is represented by a tensor product of the single-site marginals. For 
the summation vector we have the tensor representation (s\ = (1, 1)® L . 



2.2 Grandcanonical conditioning on an atypical current 

The generator H defined above is for the hopping dynamics only, it does not include 
the evolution of the integrated current J if). Since J{t) can take any integer value J, 
the state space for the full process is {0, 1} L x Z. In order to construct the matrix 
representation we choose as basis for this space the set of product vectors 1 77 ) <g) | J ) 
with JeZ. (Notice that the symbols ( • | and | • ) can be vectors in different vector 
spaces, as will be clear from the form the argument.) Defining (infinite-dimensional) 
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operators A ± through the relation A ± | J } = | J ± 1 } we obtain following Ref . [6] 
the generator 



M 

G := — ^2 [ c A cr t cr i~+i ® A+ ~ n i v i+i) + c t {a~af +1 <g> A~ - ViTii+ij] (21) 

i=-M+l 

for the full process. One sees that the elementary hopping matrices not only change 
the configuration r\ according to which jump has occurred, but also change the value 
of total integrated current J accordingly. An initial configuration is represented by 
a tensor vector \r),J) : = |?y) <S> | J). We shall take J(0) = as initial value of 
J(t). Hence an initial measure for the full process has the form | /i(0) ) (g) |0). 
In order to compute expectation values we define the extended summation vector 
( s | := ( s | (g) ( s' | where ( s' | is the infinite-dimensional summation vector for the 
counting variable J. 

Consider now the process conditioned on J(t) having reached some fixed value 
J at time t > 0. The conditional distribution /i(?7, J; t) of the ASEP under this 
conditioning has the form 

where 

Zj{t) := Prob[ J(t) = J] = (s, J|e- G '|/i(0),0) (23) 

is the marginal distribution of the integrated current and ( s, J \ := ( s \ ® ( J | with 
J G Z. We write the conditional expection of some function / of the occupation 
numbers of the ASEP as 



(f; J(t) = J) = {s,J\fe- Gt \fi (0),0) 
Zj(t) Zj(t) 



^ _ w ; °\ b ) - ° I _ \J C im< u ; ^4) 



In actual fact, however, we are not interested in conditioning on a fixed value 
J of the current, but rather in expectations for an ensemble where the current is 
allowed to fluctuate around some generally atypical value. To this end we follow 
standard procedure and define a generalized fugacity y = e s and a "grand canonical" 
current ensemble 

Y s (t) = (y^) = J2y J Zj(t). (25) 

Jez 

Analogously we define the fluctuating conditional probability of a configuration r) as 

n.(m t) = YJt) (26) 

and the corresponding expectations 

, fms _ WW®) _ Y.j&v J {f\m = j) (07 , 

uw) s - {yJ(t)) - Ya{t) ■ W 



S 



Notice that Y (t) = 1 and that therefore \ (Jq) = \ fi) is the usual (unconditioned) 
measure and ( f(t) ) is the usual (unconditioned) expectation of /. We shall refer 
to the latter quantity as the "typical" expectation. For s^Owe call corresponding 
quantities atypical. 

Consider now the unnormalized fluctuating conditional probability 

Urn t) = J2 y J( -^ J l e_G< l M°)> o )• (28) 

The following simple result is often used in the literature without proof, see e.g. [7] 
and, for a formal derivation using counting operators A ± , see [6]. One has for the 
associated unnormalized probability vector 

|/i.(t))=e-*M*|M0)) (29) 

and 

Y s (t) = (s\e-^ t \^0)) (30) 

for the normalization. Here H(s) is the matrix of dimension 2 L obtained from G by 
substituting the counting operators A^ 1 by the c-numbers e ±s . Notice that H(0) = H 
is the generator of the ASEP constructed above. For s ^ the matrix H(s) does not 
conserve probability, but nevertheless has an intuitive probabilistic interpretation. 
It gives a weight e ±s to each transition in a particular realization of the process. 
Hence we shall refer to H(s) as the weighted generator. The generalized chemical 
potential s parametrizes the mean current of this weighted process. Thus we arrive 
at the central object of our interest, which is the conditional distribution 



»s(t) ) := 



e 


-H(s)t 


Mo)) 


(« 


e -H(s)t 


Mo)) 



(31) 



This quantity describes the approach to the long-time large deviation regime from 
a given initial distribution and hence provides information about the space-time 
structure of the long-time large deviation regime. Notice that below we shall drop 
the subscript s which indicates the s-dependence. (Grandcanonical) conditioned 
expectations at time t are then computed as follows: 



(/(*)>.= 



( S \f e -H(s)t 


IMo)) 


(s 


e -H(s)t\ 


Mo)) 



(32) 



In the same spirit one can investigate the space time structure of the process 
under the condition that the integrated current across some fixed bond (k,k + 1) 
has attained a certain value. Going to the grandcanonical conditioning leads to a 
matrix H^ k \s) where only the hopping terms for the bond (k, k + 1) have the weights 
e ±s . As discussed in the introduction we refer to this setting as local conditioning, 
as opposed to the global conditioning involving the global current. Notice that 
the global time-integrated current J(t) is (trivially) extensive in system size L and 
hence its generating function does not have a good limit for L — > oo. Below we shall 
choose under global conditioning the coefficient s to be order 1/L as this yields the 
generating function for the current density j. This quantity has finite expectation 
( j ) even for L — > oo. 
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3 Shock measures 



We shall consider the evolution of two distinct types of shock measures. First we 
recall the definition of a shock measure with shock at site m in the infinite integer 
lattice Z as a Bernoulli product measure with marginal fugacity z\ up to site m — 1 
and fugacity z 2 > Z\ from site m onwards [12]. We denote these measures by 
Likewise one can define an antishock measure p~ for the infinite lattice with 
antishock at site n, i.e., fugacities z(k) = z 2 for k < n and z(k) = z\ for k > n. 

Type I shock measures for the torus are defined in analogy to these shock 
measures as follows: 

Definition 3.1 For the torus a shock measure p ! mn of type I with n ^ m is a 
Bernoulli product measure with fugacities 



z 2 at the set of sites PJ™^ 1 



{k G Tx : m < k < n} 



Zl at the set of sites P£< n :=T L \ P h 
and 



Z \ at the set of sites Pi 



{k G Tl : n < k < m} 



z 2 at the set of sites P£<™ :=T L \P l 



if-M + l<m<n<M 



if -M+l <n <m<M. 



Due to Condition S and the convention c r > q we have z 2 > Z\ and therefore 
p 2 > p\. Hence in the definition (13.1 1) the first index marks the site after which the 
high density region Pu g h begins, which extends up to site n, counted modulo L in 
the principal domain {— M + 1, . . . , M}. We call site m the microscopic position of 
the shock and site n the microscopic position of the antishock in the shock measure 
of type I. 



In vector representation we have 

®(m+M) 



with A± 



1 
1 

Z 2 



(g)(n+M) 



+ 



,2M+m-n| 



l + z 2 ] 



1 
1 

Zl 



gi(n— m) 
)(m—n) 

( 

and A~ „ 



1 

Zl 

1 

~2 



g)(iW-n) 

m < n 
i{M-m) (33) 
n < m 



[1 + z 2 ) 2M - m+n (l + Zi) 



Tensor products with exponent are defined to be absent. We remark that in terms 
of densities pi — Zi/(1 + Zi) these vectors read 



I -Pi 

Pi 
1 - P2 

P2 



®(m+M) 
(g>(n+M) 



1 - P2 

P2 
I-Pl 
Pi 



®(n— m) 



5(m— n) 



I-Pl 

Pi 
1 - P2 

P2 



5(M-m) 



m < n 



n < m. 



(34) 

The unnormalized (!) restriction of pl mn to the sector with N particles is denoted 
by ^ml i-e., 



r"m,n 



) - Pn\ pLn) 



(35) 
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where the projector on configurations with N particles is defined by 



p w |„> = ( li> » = * . (36) 

|^ (J otherwise 

Shock measures of type I are a particular microscopic realization of a shock/ antishock 
pair in a macroscopic step function density profile with density p2 in the interval 
[x,y) of rescaled coordinates m — > x, n — >■ y (modulo 1) under suitable rescaling of 
space. For M — >• oo, n — )■ oo and m fixed we recover the Bernoulli shock measures 
/z+ of [12] for the ASEP defined on Z fugacities z(k) = z\ for /c < m and ^(/c) = Zi 
for k > m. In a similar fashion one can recover an antishock measure fi~ by taking 
the thermodynamical limit such that n remains fixed and both M and m are taken 
to infinity. 

For n = M we shall drop the second subscript and write pl m := /i^M an d 
similarly for the vectors and the projections on N particles. For n = M and m = 
±M the shock measures reduce to the usual Bernoulli product measures which we 
denote by 

where y — p/(l — p) G [0, oo) is the fugacity. In particular, 

I ^m,m ) = I Mat ) = I )> I M-Af,Ar ) = I M-m ) = I )• (38) 

We stress that p} M ^ P-m- 

We shall make use of the transformation property 

Lemma 3.2 Let \ fi y ) be the vector representation of the Bernoulli product measure 
with fugacity y for L sites and N m = Ylk=m+i ^ be the partial number operator. 
Then toi, z 2 G (0, oo) and -M <m< M 

M = {t£)"~ ^">- (39) 

and for fixed particle number N 

ii \ L / \ N 

LN \ I 1 + z l \ / z 2 \ i I.N 



Proof: With the matrix representation of the number operator n for a single site 

l/*=l + (l/-l)A=(j ° 

the tensor property y Nm = l® M+m ® ^n^®M-m an( j vec ^ or representation ([33]) 
the proof of the first equality becomes elementary multilinear algebra. The second 
equality follows from N-m = N and the fact that the projection on N sites can be 
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interchanged with the number operator and that the projected vector is an eigenstate 
with eigenvalue N of the number operator N. □ 

Notice that when we work with condition S, the family of shock measures of type 
I is a one-parameter family of measures indexed by Z\. In particular, (139|) may be 
written 

(1 , \ M—m 

iT^k) q2 ^^ (41) 

and for m = —M we have the analogue of (|40l) for the sector of fixed N: 

1 + z l \ 2N\ N 



^--> = Itt^J 9 '<>• (42) 

Next we introduce a new one-parameter family of shock measures which - to our 
knowledge - has not yet been considered in the literature. 

Definition 3.3 For —M + 1 < m < M — la shock measure fi 1 ^ of type II on 
Tx with shock at position m is a Bernoulli product measure with space- dependent 
fugacities 



Ziq l for m < k < M 

Furthermore, = fij^ is a Bernoulli product measure with space- dependent fu- 
gacities 

— 2k 

z(k) = 2ig~for - M < k < M. (44) 
We write down the vector presentation of these measures 

I = ( 4-M + 1) ) 8 ( zl-M + 2) ) » ■ ■ ' ® ( ) < 45 > 

with 

m M 

Z m = J] (l + V 2 ^) II (l + ^i? 212 ^) (m<n). (46) 

k=-M+l k=m+l 

The unnormalized restrictions of these measures on the sectors with N particles are 
denoted 

Shock measures of type I and type II are related by a similarity transformation. 
This is the content of the following lemma. 

Lemma 3.4 Assume condition S and let 

Then for —M < m < M we have 

l/40 = ^^%K> (48) 
12 



with 



Y 



n 



1 + zi<r L 

1 + z x 



k=-M+l 

Moreover, for fixed particle number N 



u n > N 



1 

y~ 

1 n 



M 

n 

k=m J r l 



o m-f- M — k 

1 + z x q L 
1 + g 2 zi 



2(m-M) . r 



(49) 



(50) 



Proof: Observing that z^q 2 ~ = z^q 2 = z 2 the proof works like for lemma 
but using also that not only the projector P^ on N particles and the number operator 
N are diagonal in the basis used throughout this work, but also the transformation 
U. Hence all three operations can be arbitrarily interchanged. □ 

To elucidate the nature of these type II shock measures we define parameters c, 
E through e 2E := q 2 , e~ 2Ec := Z\. Using Condition S gives z 2 = e~ 2E ^ c ~^ and we 
can then express n 1 ! as a Bernoulli product measure with densities 



p{k) 



tanh 



E(k+cL) 



1 - tanh f 



—M < k <m 
m < k < M 



(51) 



For hopping asymmetry q of order 1 the variation of density between neighbouring 
sites is of order 1/L, except between site m and m + 1 where there is a jump of order 
1. On the other hand, for strong asymmetry with E of order L the density jumps 
from nearly to nearly 1 at the shock position m and there is a sharp transition from 
nearly 1 to nearly 0, extending over a few lattice sites in the vicinity of k* = —cL 
(modulo L). Here "nearly" means up to corrections exponentially small in L. This 
sharp transition is the microscopic analogue of an antishock in the type II shock 
measure which exists only for strong asymmetry. For small asymmetry the density 
profile is almost linear, see Figure [TJ 



4 Microscopic conditioned evolution of shocks 
4.1 Main results 

We are now in a position to state and prove our main results on the time evolution 
of the ASEP under conditioned dynamics, viz. for type II shocks under global 
conditioning, and for type I shocks under local conditioning. Throughout this section 
we assume condition S to be satisfied. 

Theorem 4.1 (Shock evolution under global conditioning) Consider the ASEP on 
Tl with N particles and hopping rates q and c r to the left and to the right respec- 
tively. Let the initial distribution at time t = be given by the type II shock measure 
(j>m' N and let ^^(t) denote the grandcanonically conditioned distribution of the 
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Figure 1: Interpolated density profiles for type II shock measures for lattice size 
L=200, shock position m — 10 and different asymmetries and total densities: (a) 
E = 0.2, c = 0.5, (b) E = 2, c = 1, (c) E = 2, c = 0.5, (d) E = 40, c = 0.3. 
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ASEP at time t under global conditioning with weight e sL = Ci/c r . Then for m G T L 
and any t > 

M 

/#"(*) = E Vt{l\m)^ N (52) 

i=-M+l 

where 

°° / r \ ( l ~ m +P L )/ 2 
p t (l\m) = e-^ + ^ t £ (^J I m ^ l+pL (2y^J £ t) 
p=— 00 ^ ' 

wit/i i/ie modified Bessel function I n (-) and c r = c r q~^~ and q = c^q^~ . 

We remark that p t (l \ m) is the probability that a particle that performs a 
continuous-time simple random walk on with hopping rates c r to the right and 
q to the left, respectively, is at site I at time t, starting from site m. Hence the 
interpretation of (152]) is that at any time t > the distribution of this process is a 
convex combination of shock measures of a conserved microscopic structure equal to 
that of the initial one. The probabilistic interpretation of the weights in this convex 
combinations allows us to say that the shocks (or better, the shock positions) per- 
form continuous time random walks on with hopping rates c r (to the right) and 
q (to the left). Seen from the shock position the measure is invariant. As opposed 
to the macroscopic results of [8] our results provide exact information for any t > 
and for finite lattice distances. Notice however, that our results require the choice 
e sL = ce/c r with arbitrary asymmetry q, while the phenomena observed in [S] are 
valid for any s, but require weak asymmetry. 

Proof of Theorem In order to prove (1321 we observe that with (T4TT) (which 

follows from (1381) and (1391) in Lemma (I3.2p ) and with Lemma (13 .4p we can write the 
initial measure 

\^' N ) = V m \^) (53) 

with 

V m = —qx(^T=-M+ii-M+m-k)n k +j:^ =m+1 (M^-k)h k ) ^ 
and a non-zero normalization factor 

M 

W = (l + Zl ) L H (l + z iq - 2 ^ (55) 

k=-M+l 

that does not depend on m. We recall ([H]) and get for the unnormalized measure 

^ t \^: N (t)} = -v m H^\^ 1 (t)) (56) 

with = V- x EV m . 

Similarity transformations of the diagonal number operators on the particle cre- 
ation and annihilation operators are of the form 

p -a k h k ± a k h k _ j °\ k I , 

aiG -{ e^o-f k = l ■ (57) 
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From (1571) and Condition S we thus find 

M 

ff(m) = _ [c r {afa~ +1 - hiV i+1 ) + c e (a^cr+ +1 - ^+1)] (58) 

i=-M+l 

~ [ c l a m a m+l ~ c rh m Vm+l + Cr^^m+l ~ CiVrnfim+l] (59) 

where the prime at the summation indicates the absence of the term with i — m. 
We add — (c r — cg)(n i+ i — hi) to each term in (155]) and — (c r — ce)(h m+ i — h m ) to 
(|59p . Since the sum of these terms is zero we arrive at 

M 

ffM = _ [c r {afa- +1 - Vih i+1 ) + c £ {a-af +1 - hiV i+1 )} (60) 

i=-M+l 

- [ci(a+a~ +1 - h m v m+1 ) + c r (a~a+ +1 - v m h m+1 )] (61) 

which we write in the form if ( m ) = + f?( m ) with bulk term iJ^ = ^ h\ given 
by (16"U|) transformation term _E?( m ) given by (1611) . 

Using the properties of the particle creation and annihilation operators listed 
above it is easy to check by straightforward computation that h\\ /i^ ) = 0. On the 
other hand, B^\^ > = -(c r - c e )(h m+1 -h m )\^). This implies V m B^\ ^ ) = 
— (c r — ci)(h m+ i — h m )V m \ } since V m and the number operators hi are both 
diagonal and hence commute. With the projector property h\ = h^ one has 
q- 2N ' L V m+1 = (1 + (q~ 2 - l)h m+1 )V m and q 2N l L V m . x = (1 + (q 2 - l)h m )V m . There- 
fore c r q~ 2N/L V m+1 + Qg 2Ar/L ^ m _i - (c r + cg)V m = -(c r - c £ )(h m+1 - h m )V m . Putting 
these results together yields 

^ N (t) = crq-^i&dt) + W™tt- N i(t) ~ (c r + c e )tt> N (t) (62) 

for any m ET L and any t > 0. 

Next we compute the time derivative of the normalization R m (t) = ( s \e~ Ht \ )■ 
Using (1621) we get 

±R m (t) = -(s\e m H\^ N ) 

= c r q- 2N ' L R m+1 (t) + c e q 2N / L R m ^(t) - (c r + c,)R m {t) (63) 
= [c r q- 2N ' L + m 2N ' L - ( Cr + Ci)] Rm(t) 

where the last line follows from translation invariance. Therefore 

Vk' N (t) = ti I,N (t)/R(t) = exp [(-c r q- 2N / L - Cl q 2NIL + c r + c^t}^ (t) / R(0) 

(64) 

which implies the system of linear ODE's 

^(t) = C r q- 2 ^ L ^(t) + C^^it) - {c r q- 2N ' L + Cl q 2NIL ) ^ N (t). 
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Now it only remains to show that (152)) satisfies this system of ODE's with the initial 
condition /x„(0) = pJ^. It is elementary that the random walk transition probability 
of the theorem satisfies the forward evolution equation 

±p t (l | m) = c r q- 2N ' L p t {l-l | m ) +Ce q 2N / L p t (l+l \ m)-{c r q~ 2N l L + Cl q 2N l L ) p t (l \ m) 

(66) 

with initial condition p (l \ m) = 8^ m . The theorem thus follows from translation 
invariance p t (l + r | m + r) = p t (l \ m) Vr G Z and periodicity p t (l+pL \m) = p t (l\ m) 
Vp G Z of the transition probability. □ 

Theorem 4.2 (Shock/ Antishock evolution under local conditioning) Let the initial 
distribution of the ASEP with N particles be given by the type I shock measure /i^ 
satisfying condition S with shock at bond m and antishock at bond M and let p^ N (t) 
denote the unnormalized grandcanonically conditioned measure at time t under local 
conditioning at bond M with weight e s = ci/c r . Then we have V£ > 

M 

ffi(t)=^-*- e # Vt{l\m)8 l - m q m ^ d N (67) 

l=-M+l 

with p t {l | m) defined in (fJT7]j and 5 = 1+ + 2^ 1 ■ 

Proof: First we consider the evolution of unnormalized measures jx 1 ^ (t). By con- 
struction (described informally in Sec. 2) the matrix H defined by (158]) and bound- 
ary term Bm given by (1591) is the generator for the grandcanonically conditioned 
evolution under local conditioning at bond (M,—M + 1) with weight e s = c £ /c r . 
Therefore d/dt| fa N ) = -H M \ jiff). Notice next that U = WV M where U is de- 
fined in (I47p and Vm is defined in f[5~4"l) . Therefore we may write H M = U~ l HU and 
we have by Theorem (14.11) in conjunction with Lemma (13.41) 



uM\~I,N\ _ v n'^^^T^^TT- 1 U\ T,U,N \(+\ 



Y a 2 



+c £ q~ 



'-N 



c r q 



_ 2 M-rn-l N J N 



Y 



m+l 



_ 2 M-rn + l N „J N 



1 m—1 



{c r + c e )—q ' 

* m 



2 M=ui N , ~ j N 



Am ) 



f'm+l I 



-M <m< M). (68) 



From (j4"9~|) we have Y m /Y m+ i = (1 + Zi)/(1 + z 2 ) and therefore 

d ti?{t) = c r 5~S m N +l (t) + ceS-^it) ~ (cr + cdtffQt) Vm g T \ {—M + 1, M}. 

(69) 



dt 
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For m = M we have 

rM\~I,N\/.\ xr rr-l £rl r.H,N 



-H M \~^){t) = -Y M U^H\~^)(t) 
= Ym 



-M 1 _ 2 ^i TV I ~I,N \ 



Af+l 



. 2iv 1 _ 2 J-/Vi ~I,N \ 

J- m 



(70) 



Straightforward computation yields Y M /Y_ M+1 = [(1 + z 2 )/(l + <2i)] L 1 - 
For m = — M we use Lemma (13. 2p . This yields 

^M^*) = C r^ L Q~ 2N ^-AI+l( t ) + C^ _1 /iM-l(0 - (Cr + Q^lfW (71) 

^A-Af+iW = Cr^-l +2 (*) + c e 5~ 1 5 L q 2N fi I ^(t) - (c r + q)/^M +1 (i). (72) 

From (155|1 - (1721) one finds that he transformed measure | ) (t) = S m q 2Nm / L | ) (£) 
satisfies the system of linear ODE's (16"2l . Hence, by the same arguments as used in 
flUD we obtain ([67]). □ 

Corollary 4.3 T/ie normalized grandcanonically conditioned measure at time t of 
Theorem \4-W A 4 ™^) ^ s of the ft 



orm 



M 

W®=C® E ftGI™)* 1 "*^^ (73) 



Z=_M+1 



wii/i a normalization constant C(t) = 'YldL-M+iPt{l\ m )b~ l rn Q <L '(^l/ 1 ^) that 
is strictly positive and finite for allt > 0. In particular, 



Z=-A/+l 

ttfif/i C* = Xlz^-Af+i ^ l ~ m q x ~ ( s I Z-4 ,7V ) ^ e unique stationary conditioned lim- 
iting measure. 

The properties of C(t) are obvious from the fact it is a transition probability for 
a random walk with unnormalized positive initial distribution and that the random 
walk propagator on satisfies lim^oo Pt(l \ m ) = l/LVl,m G Tl- The unique- 
ness of the measure follows from the Perron- Frobenius theorem for the weighted 
generator. □ 
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5 Final remarks 



Our results are exact results on the lattice scale and for finite times. With regard 
to the macroscopic large- deviation theory of [8] Theorem (14.11) suggests that there 
exists a value of s(v) such that the optimal macroscopic density profile obtained 
from the large- deviation theory of [5] in the limit v — > oo is of the form p(x, t) = 
1/2[1 — tanh (E(x + a — vt))} with speed v = c r exp (— 2ps) — Qexp (2ps) and a jump 
discontinuity at some point x*(t) that also travel with speed v. The stationary 
solution in Theorem f!4.2[) for local conditioning with weight e s = cijc r indicates the 
existence of an antishock that remains microscopically sharp and does not fluctuate. 
This behaviour is in sharp contrast to the typical behavior of the ASEP where an 
antishock evolves into a rarefaction wave. 

The evolution of the shock initial measures as a one-particle random walk indi- 
cates a self-duality relation of the conditioned ASEP with periodic boundary con- 
ditions, which is absent for the unconditioned dynamics. Results from [TJ] for the 
microscopic dynamics of more than one shock, each satisfying condition S, and prop- 
erties of the underlying representation theory of the quantum group U q [SU(2)} for 
periodic systems [21] suggests that one may study the behaviour of shock/ antishock 
measures with several shocks under global and local conditioning. One then expects 
the evolution of n shocks to be given by the dynamics of n interacting particles. 
The results of [13] might then provide information about more general solutions of 
the large deviation theory and about their fluctuations. 
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